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Fermionic Molecular Dynamics (FMD) models a system of fermions by 
OO ' means of many-body states which are composed of antisymmetrized prod- 

ON ' ucts of single-particle states. These consist of one or several Gaussians 

on : localized in coordinate and momentum space. The parameters specifying 

them are the dynamical variables of the model. As the repulsive core of the 
C"| \ nucleon-nucleon interaction induces short range correlations which cannot 

. be accommodated by a Slater determinant, a novel approach, the unitary 

correlation operator method (UCOM), is applied. The unitary correlator 
OO , moves two particles away from each other whenever their relative distance 

is within the repulsive core. The time-dependent variational principle yields 
the equations of motion for the variables. Energies of the stationary ground 
states are calculated and compared to exact many-body results for nuclei up 



> 

o 



to 48 Ca. Time-dependent solutions are shown for collisions between nuclei. 



O ! 1 Concept 

OO n 

ON ■ The general concept of Fermionic Molecular Dynamics (FMD) H is to describe a system of 

fermions by a many-body trial state which is antisymmetric with respect to particle exchange 
("Fermionic"). The attribute "Molecular" indicates that the individual single-particle states, 
from which Slater determinants are formed, are wave packets localized in coordinate as well 
as momentum space. Hence they are the quantum analogue to classical points in phase space. 
The time evolution of the parameters which specify the many-body state, like mean positions, 
mean momenta, or spin orientations of the wave packets, is obtained from the time-dependent 
variational principle ("Dynamics"). A closely related model is AMd! which follows the same 
concept but differs from FMD in several aspects. 

1.1 The trial state 

In FMD the trial states are antisymmetrized A-ho&y states 

\Q) = CA \qi) ® \q 2 ) ® • • • ® \u) , (1) 

where Q denotes the set of all subsets tjk which specify the parameters of the single-particle states 
\q~k)- These are Gaussians, localized in phase space, or superpositions of several Gaussians. In 
this case q^ includes in addition to the complex parameters mean positions, widths and two- 
component spinors also the complex amplitude for each Gaussian. In general \Q) can be a 
superposition of correlated Slater determinants, then the set Q includes also the configuration- 
mixing amplitudes. 

A is the projector onto the antisymmetric Hilbert space of fermions. C is a unitary correla- 
tion operator which takes care of the short range correlations as explained in subsection 0. 

a Invited talk at "XVII RCNP International Symposium on Innovative Computational Methods in Nuclear Many- 
Body Problems", November 10 - 15, 1997, Osaka. To be published in the proceedings, World Scientific. 
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1.2 The time-dependent variational principle 



The equations of motion for all parameters in Q = {q v ; v = 1, 2, • ■ • , N} are obtained from the 
time-dependent variational principle 

5 j\t(Q(t)\ i^-H \Q(t)) = 0, (2) 

where the variation is taken with respect to all individual parameters q^ or their complex con- 
jugate q* The resulting equtions of motion are 

_ B !««)> = (^OWlE*^ _ H m) = o. (3) 

The physical meaning of these equations is that the deviation of the exact state from the ap- 
proximate \Q(t + dt)), which develops during a time interval dt, 

\Error(t + dt)) = (1-idtH) \Q(t)) - \Q(t + dt)) (4) 

= [dt(iY,^w- - H )\Q(t)) (5) 

is orthogonal to the tangent Q(t)\ ^ = Q(t)) = g|- \Q(t)) at each time t, as can be 
seen by comparing Eq. (||) with (|3|). In other words the variational principle chooses the best 
\Q(t + dt)) within the manifold \Q). 

Stationary solutions of the variational principle (^), ground states for example, imply 

H \Q(t)) = 0. (6) 

dq* 

1.3 Quantum-branching 

If one wants to allow the time evolution to deviate from \Q(t)) without increasing the degrees 
of freedom one may take the operator 

AV(t) = ( ij>5£- — H ) (7) 

from Eq. @ as a time-dependent pertubation to calculate transition amplitudes for quantum- 
branching away from \Q{t)). If the time evolution within the manifold happens to be exact, 
AV(t) \Q(t)) = and no branching occurs because all quantum physics is already in \Q(t)). 
The branching introduced by Ono and Horiuchi in AMD-VeI or the branching between energy 
eigenstates proposed by Onishi and RandrupB may be investigated under this aspect. 



2 The Unitary Correlation Operator Method (UCOM) 

In nuclear physics this straightforward concept encounters a complicated Hamiltonian. Due to 
the complex structure of the nucleons and the mesons exchanged between them, the interaction 
depends on spin and isospin, is momemtum dependent and overall repulsive at short distances, 
see for example Ref. H. The properties of the Hamiltonian are of course intimately related to 
the choice of the trial state. The tensor part induces a strong correlation between the spins of a 
pair of nucleons and the direction of their relative distanced. The repulsion at short distances 
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reduces the probability amplitude to find two nucleons close together. Both correlations cannot 
be described by a (antisymmetrized) product state. Therefore the most simple ansatz of a single 
Slater determinant has to be modified in order to incorporate the above mentioned correlations. 
In the following we sketch the new Unitary Correlation Operator Method (UCOM) which allows 
an approximate treatment of the short range repulsion even in time-dependent states. 

Due to the short-ranged repulsive core in the nucleon-nucleon potential V the many-body 
state is depleted as a function of the relative distance | for each pair (ij) when they 

are close to each other. These short range correlations cannot be described by shell model states. 
The most common procedure to remedy this problem is Brueckner's G-matrix method which 
replaces the bare V by an effective interaction G. Another method is the Jastrow approach 
where the correlated ground state of the nucleus is assumed to be of the form Y\i<j f( x ij) 1^) 
where f(xij) is a correlation factor which diminishes the probability to find two nucleons at 
small distances x\y We propose a new method i in which the correlated state \Q) is obtained 
by a unitary transformation C = exp{— iG}. With respect to the choice of G it differs from 
the work of K. Suzuki et all, where an effective interaction for the shell model is defined. In 
collisions between nuclei a preferred shell model space does not exist. Therefore, our definition 
of G is not related to a given basis but to the behaviour of the relative wave function at small 
distances. The correlated state is given by 

\Q) = C\Q)=e- [G \Q) , (8) 



where the generator G is a hermitean two-body operator which depends in general on the relative 
distance Xij, the relative momentum q iv the spins and isospins of the two nucleons 0. For sake 
of simplicity we sketch here only the case of a repulsive core without spin-isospin-dependence. 
In the application shown later G depends on spin and isospin such that the different channels 
of a central interaction can be correlated individually. 

The aim of C is to push two nucleons away from each other whenever they get too close. 
The most simple ansatz which does that is 

G = -k Yl \ ( <iiizr- I s 0%) + s ( x ij) I —'hi ) \ , = \xij\ . (9) 





s(xij) is roughly speaking the distance which the particles i and j are moved away from each 
other by exp{— iG} if they are found at a distance x%j in \Q). s(xij) is largest if Xij lies inside 
the hard core and s(xij) — ► if x^ is outside the repulsive interaction. 

As = C _1 is unitary one can correlate the states or equivalently the observables: 

(Q\H\Q) = (Q\C^HC\Q) = {Q\C- l HC\Q) = (Q\ H \Q) . (10) 

Closed expressions for the correlated relative wave function, the correlated potential and the 
correlated kinetic energy can be given by using the coordinate transformations Xij — > R±(xij) 
defined by 

rR±(x) (j£ 



.(0 M - 

For example the correlated potential is simply 

vW( Xij ) = C- l V( Xij ) C = V(C- 1 x lj C) = V(R+{ Xij )) . (12) 
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Figure 1: Uncorrelated state (x\2\4>) (dashed line), correlated state (S12 10) (full line), Afnan-Tang potential 
V(|cci2|) (thin line) and correlated potential y' 2 '(|Si2|) (thin dashed line) as a function of the relative distance 
X12 = |2ia|. 

Or a correlated relative wave-function can be expressed in terms of R- which is the inverse of 

R_^_ clS i^X - — - X ^ *^ j — "^ij ) 

(x\cf>) = (x\C\$) = ^MJ±R_( X ) (R-(x)-\$) (13) 
x V ax x 

Fig. 1 displays the radial dependence of the correlated and uncorrelated deuteron wave 
function together with the Afnan Tang S3 potentials. The functional form of s(x) or equvalently 
of R+{x) is chosen such that the correlated state \4>) for short distances equals the exact solution. 
The correlated kinetic energy is more involved as one has to calculate the correlated relative 
momentum C -1 ^ C. The details are given in Ref.lZI. 

Once the correlator C is adjusted to reproduce the two-body system at low energies (long 
wave length) one can calculate the ground state energies of nuclei by minimizing (Q\H \Q) = 
(Q \ H\Q) with respect to \Q) which in this application is a single Slater determinant composed 
of localized Gaussians (FMD). 

Unlike the original Hamiltonian H = T + V, the correlated Hamiltonian 

H = C~ X H C = C~ X T C + C~ X V C = T +f ^ +V [2] +f [3] +V [3] + ■■■ (14) 

is not a one- plus two-body operator anymore. It contains three-body and higher interactions 
because the generator G is a two-body operator. 

We calculate and V®, which are functionals of R+, analytically (for example as in 
Eq. (|l^)) and neglect irreducible three-body and higher terms in the expansion (|T4|). This turns 
out to be a very good approximation at typical nuclear densities. Estimations of the three-body 
terms give corrections less than 5% of binding energy for the a- particlel. 

Fig. 2 compares uncorrelated (left columns) and correlated (right columns) energies. The 
correlated potential energy (Q|l/t 2 l \Q) (grey bars at negative values) is about twice the un- 
correlated (Q\V \Q) in all nuclei. This gain in binding is counteracted by an increase in the 
kinetic energies (dark grey bars at positive values). Both together yield binding energies (black 
bars) which are within 8% deviation from results of Yakubovski calculationslll for 4 He and CBF 
calculations^! for "O and <^Ca. This sunrises since H C2 = T+f^ + V^ is the same for all 
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Figure 2: Energies for the modified Afnan-Tang S3 potential. Left columns: expectation values of kinetic energy 
{ T ) , potential { V } and sum of both { H ) . Right columns: expectation values of correlated kinetic energy 
( T + f [2] } , correlated potential { V^ 21 } and sum of both { H C2 ) . Arrows indicate the results of a Yakubovski 
calculation for 4 He and of CBF calculations for 16 and for 40,48 Ca. 



nuclei and not density dependent but it contains momentum dependent parts in the two-body 
part of the correlated kinetic energy . 



3 Ground states with several Gaussians 

The Gaussian single-particle states do not describe properly the surface of nuclei because their 
tails do not fall off exponentially. Especially in the case of halo nuclei this is an important part 
of the physics. As Gaussians represent an overcomplete basis, any shape can be achieved by 
superimposing several Gaussians. In the following example we minimize the energy of 6 He for 
a single Slater determinant with one Gaussian per particle and with two co-centered Gaussians 
with free widths and relative amplitudes. The resulting proton and neutron densities are dis- 
played in Fig. 3 on a logarithmic scale. For the one Gaussian case (l.h.s.) the distribution of the 
neutrons is, as expected, wider than for the protons, but already with two Gaussians (r.h.s.) the 
trial state is much improved and exhibits nice exponential tails of different ranges. The binding 
energy increases by about 3.5 MeV and the 4 He core which is distorted on the l.h.s. assumes 
again equal proton and neutron densities in the center. Three Gaussians do not improve the 
picture anymore. This example shows that a Gaussian basis may be a good representation for 
nuclear structure calculations. See for example also the contribution by M. Kamimura in these 
proceedings. 

We have also diagonalized the Hamilton matrix (Qk \ H \Qi) within a subspace spanned by 
the nonorthogonal many-body states \Qk)- Dimensions of up to several hundred are feasible 
and for example eigenstates of total spin may be obtained by taking \Qk) as rotated intrinsic 
states. See also contributions of S. Aoyama, N. Itagaki and S. Okabe or Y. Kanada-En'yo and 
H. Horiuchi or of K. Varga and Y. Suzuki in these proceedings. 
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Figure 3: Proton and neutron densities of 6 He with one Gaussian per aprticle (l.h.s.) and two Gaussians per 
particle (r.h.s) 



4 The time dependent case 

In the time-dependent case one cannot be so ambitious as for ground states because all deriva- 
tives in Eq. |3| have to be calculated at each time step. Here we have to compromise and take 
only one Gaussian per particle and one Slater determinant. In addition the Hamiltonian is a 
simplified version of H = C~ 1 HC. We tried to adjust the parameters of the central H such 
that ground state energies of nuclei up to mass 56 are reproduced within about 5% error. In 
doing that we found several Hamiltonians which describe similarly well the ground states but 
differ in their momentum dependence. A particular H has two almost degenerate local energy 
minima in the manifold \Q). One minimum has states with an intrinsic cluster structure with 
narrow Gaussians placed at different positions in coordinate space, but centered in momentum 
space. The other shows co-centered Gaussians in cordinate space which are slightly displaced in 
momentum space, which resembles a shell model state. In Fig. || one sees that the two different 
types of states behave quite differently when they collide. In the left column, where two cluster 
states interact, multifragmentation is observed which could be characterized as shattering of the 
nuclei. Shell model states on cluster states produce also several fragments but in addition many 
nucleons are emitted. They are seen in the middle column as crosses (centroids of the packets) 
without surrounding contour lines, because the widths are already so large that the density 
distribution of the emitted nucleon has fallen below the lowest contour line. Finally, reactions 
with two shell model states (right column) look more like a deep inelastic collision with many 
nucleons evaporated because of the high excitation energy. 
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Figure 4: Densities projected on scattering plane for Ca+ Ca collisions at b=2.75 fm and Ei a b=S5 AMeV. 
Initial nuclei in cluster and shell-model configurations. Crosses denote centroids of wave packets. 
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